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Abstract. We study the propagation, observation and control properties of the 1 — d wave equation on 
a bounded interval discretized in space using the quadratic P2 - classical finite element approximation. A 
careful Fourier analysis of the discrete wave dynamics reveals two different branches in the spectrum: the 
acoustic one, of physical nature, and the optic one, related to the perturbations that this second-order finite 
element approximation introduces with respect to the Pi one. On both modes there are high frequencies 
with vanishing group velocity as the mesh size tends to zero. This shows that the classical property of 
continuous waves of being observable from the boundary fails to be uniform for this discretization scheme. 
As a consequence of this, the controls of the discrete waves may blow-up as the mesh size tends to zero. 
To remedy these high frequency pathologies, we design filtering mechanisms based on a bi-grid algorithm 
for which one can recover the uniformity of the observability constant in a finite time and, consequently, 
the possibility to control with uniformly bounded L 2 - controls appropriate projections of the solutions. 
This also allows showing that, by relaxing the control requirement, the controls are uniformly bounded and 
converge to the continuous ones as the mesh size tends to zero. 



1. Introduction, problem formulation and main results 
Consider the 1 — d wave equation with non-homogeneous boundary conditions: 

( y tt (x, t) - y xx (x, t) = Q, x € (0,1), t > 0, 

(1) \ y(o,t) = o, y(l,t) = v(t), i>0, 

[ y(x,0) =y°(x), y t (x,0)=y 1 (x), x G (0, 1). 

System (JlJ is said to be exactly controllable in time T > 2 if, for all (y°, y 1 ) eL 2 x /f _1 (0, 1), there exists a 
control function v £ L 2 (0,T) such that the solution of Q can be driven at rest at time T, i.e. 

(2) y(x,T) = y t (x,T)=0. 

We also introduce the adjoint 1 — d wave equation with homogeneous boundary conditions: 

{Uu(x, t) - u xx (x, t) = 0, x G (0, 1), t > 0, 
u(0,t) = u(l,t) = 0, t>0, 
u(x,T) = u°(x), Ut(x,T) = u^x), x G (0,1). 

This system is well known to be well posed in the energy space V := Hq x L 2 (0, 1) and the energy below 
is conserved in time: 

£{u°y) :=i(|K,t)||^ (0il) + |K(-,t)||| a{0 , 1) ). 

The Hilbert Uniqueness Method (HUM) introduced in [2] allows showing that the property of exact 
controllability for (JlJ is equivalent to the boundary observability of ^3j. The observability property of the 
wave equation ensures that the following observability inequality holds for all solutions of Q provided T > 2: 

T 

|2 



(4) £(r/V) < C(T) J\u x (l,t)f 



dt. 
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The best constant C(T) in Q is the so-called observability constant. The observability time T has to be 
larger than the characteristic time T* := 2 which is needed by any solution associated to initial data (u°, u 1 ) 
supported in a very narrow neighborhood of x = 1 to travel along the characteristic rays x{t) = x — t, get to 
the boundary x = and come back to the boundary x = 1 along the characteristics x(t) = x + 1. 

It is also well known that for all T > and all solutions u of the adjoint problem ^ with initial data 
(v?,v}) € V, the following admissibility inequality holds: 

T 

(5) c{T) [ \u x {l, t)\ 2 dt<S{u ,v}), 



o 



so that for all T > 2, Hc^^l, -)Hl 2 (o,t) and 1 1 xt 1 ) 1 1 v are equivalent norms. 

As a consequence of these results, it is easy to see that for all (y°, y 1 ) s L 2 x if _1 (0, 1), there exists a 
control v € L 2 (0,T) driving the solution of Q to the rest at t = T, i.e. such that ^ holds. This turns out 
to be equivalent to the fact that 

T 

(6) / v(t)u x (l, t) dt = (u(-, Q) ,«*(•, Q)))v,v, 



o 

for all solutions u of where (•, -)v',v is the duality product between V' = H^ 1 x L 2 (0, 1) and V. 
The HUM control v, the one of minimal L 2 (0, T)-norm, has the explicit form 

(7) v(t)=v(t):=u x (l,t), 

where u(x, t) is the solution corresponding to the minimum (u°,u 1 ) G V of the quadratic functional 

T 

(8) J(u°,u l ) := \J \u x (l,t)\ 2 dt- ((y\-y o ),(u(;0),u t (;0)))v,V- 



The effect of substituting the continuous controlled wave equation ([I]) or the corresponding adjoint 
problem ([3| by discrete models has been intensively studied during the last years, starting from some simple 
numerical schemes on uniform meshes like finite differences or linear P± - finite clement methods in [13 and, 
more recently, more complex schemes like the discontinuous Galerkin ones in [17] . In all these cases, the 
convergence of the approximation scheme in the classical sense of the numerical analysis does not suffice 
to guarantee that the sequence of discrete controls converges to the continuous ones, as one could expect. 
This is due to the fact that there are classes of initial data for the discrete adjoint problem generating high 
frequency wave packets propagating at a very low group velocity and that, consequently, are impossible to 
be observed from the boundary of the domain during a finite time, uniformly as the mesh-size parameter 
tends to zero. This leads to the divergence of the discrete observability constant as the mesh size tends to 
zero. High frequency pathological phenomena have also been observed for numerical approximation schemes 
of other models, like the linear Schrddinger equation (cf. llj), in which one is interested in the uniformity 
of the so-called dispersive estimates, which play an important role in the study of the well-posedness of the 
non- linear models. 

Several filtering techniques have been designed to face these high frequency pathologies: the Fourier 
truncation method (cf. [13] ). which simply eliminates all the high frequency Fourier components propagating 
non- uniformly; the bi-grid algorithm introduced in [5] and [5] and rigorously studied in |12j . [16) or |20) in the 
context of the finite differences semi-discretization of the 1 — d and 2 — d wave equation and of the Schrddinger 
equation (cf.|llj). which consists in taking initial data with slow oscillations obtained by linear interpolation 
from data given on a coarser grid; and the numerical viscosity method, which by adding a suitable dissipative 
mechanism damps out the spurious high frequencies f|19j. |21jV We should emphasize that the mixed finite 
element method analyzed in [5j is, as far as we know, the unique method that preserves the propagation 
and controllability properties of the continuous model uniformly in the mesh size without requiring any 
extra filtering. The interested reader is referred to the survey articles [7] and |22j for a presentation of the 
development and the state of the art in this topic. 
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The purpose of the present paper is to analyze the behavior of the quadratic Pi - finite element semi- 
discretization of the problems ([T]) and ([3]) from the uniform controllability/observability point of view. In 
Section [2] we introduce in a rigorous way the discrete analogue of ([I]) and ^ and explain the minimization 
process generating the discrete controls. In Section [3] we analyze the spectral problem associated to this 
discrete model and reveal the co-existence of two main types of Fourier modes: an acoustic one, of physical 
nature, related to the nodal components of the numerical solution, and an optic one, of spurious nature, 
related to the curvature with which the quadratic finite element method perturbs the linear approximation. 
We also study finer properties of the spectrum, for example the spectral gap, identifying three regions of 
null gap: the highest frequencies on both acoustic and optic modes and the lowest frequencies on the 
optic one. The content of this section is related to previously existing work. For instance, the dispersive 
properties of higher-order finite element methods have been analyzed in [T] in the setting of the Hclmholtz 
equation. An explicit form of the acoustic dispersion relation was obtained for approximations of arbitrary 
order. It was shown that the numerical dispersion displays three different types of behavior, depending on 
the approximation order relative to the mesh-size and the wave number. In Section [4] we obtain some 
spectral identities allowing us to analyze the discrete observability inequality for the adjoint system. In 
Section [5] we show that the Fourier truncation of the three pathological regions of the spectrum leads 
to an uniform observability inequality. In Section [6] we prove that a filtering mechanism consisting in, 
firstly, considering piecewise linear initial data and, secondly, preconditionning the nodal components by a 
bi-grid algorithm guarantees uniform observability properties. Within the proof, we use a classical dyadic 
decomposition argument (cf. |12|). which mainly relies on the fact that for this class of initial data the total 
energy can be bounded by the energy of the projection on the low frequency components of the acoustic 
dispersion relation. We should emphasize that our results are finer than the ones in [3] or [5], where one 
obtains uniform observability properties for finite element approximations of any order, but by filtering the 
Fourier modes much under the critical scale l/h. Here we only consider the particular case of quadratic 
finite clement approximation on 1 — d meshes, but we get to the critical filtering scale l/h. Note however 
that the results in [3] and [3] apply in the context of non-uniform grids as well. In Section [7] we present 
the main steps of the proof of the convergence of the discrete control problem under the assumption that 
the initial data in the corresponding adjoint problem are filtered through a Fourier truncation or a bi-grid 
algorithm. Section [8] is devoted to present the conclusions of the paper and some related open problems. 

2. The P 2 - finite element approximation of 1 — d waves 

Let N S N, h = l/(N + 1) and = xq < Xj < xn+i = 1 be the nodes of an uniform grid of the 
interval [0, 1], with Xj — jh, < j < N + 1, constituted by the subintervals Ij — (xj, Xj+x), with < j < N. 
We also define the midpoints Xj +1 / 2 = (j + l/2)h of this grid, with < j < N. Let us introduce the 
space Vp(a,b) of polynomials of order p on the interval (a, b) and the space of piecewise quadratic and 
continuous functions Uh :— {u G Hq(0, I) s.t. u\j. 6 "Pi{Ij)t < j < N}. The space Uh can be written as 
Uh = span{<^-, 1 < j < iV}©span{0 J+1 / 2 ,O < j < N}, where the two classes of basis functions are explicitly 
given below (see Fig. [I]) 

!j^(x - Xj-i/ 2 )(x - Xj-\), X £ Ij_x, 
w {x- x j+1/2 )[x - x j+ i), x£lj, and <p j+1/2 (x) = 
0, otherwise 

The quadratic approximation of the adjoint problem ^ is as follows: 

( 9 ) j Find u h {-, t) € U h s.t. ^(u h (-,t),(p) L 2 (oX) + (u h (-,t),(p) H i {0A) = 0,yip eU h , 
\ u h (x,T) = ul(x),u hjt (x,T) = u\(x), ace (0,1). 

Since Uh(-,t) £ Uh, it admits the decomposition Uh{x,t) — Y^j^i 1 Uj/2{t)<j)j/ 2 {x). The function uj-^^t) can 
be identified with the vector of its coefficients, \Jh(t) := (Uj/2(t))i<j<2N+i (in the sequel, all vectors under 
consideration will be column vectors). Thus, using ip — <f>j/2, 1 < J < 2iV+ 1, as test functions in (l9|, system 
([9]) can be written as the following system of 2N + 1 second-order linear differential equations (ODEs): 

(10) M h U h>tt (t) + S h U h {t) - 0, V h (T) = U°, U /M (T) = u£, 



^(x- Xj)(x- x j+ i), x e Ij, 



0, otherwise. 
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X H X i-1/2 X i X j+1/2 X j+1 




Figure 1. The basis functions: <j)j (left) and ^j+1/2 (right). 



where Mh and Sh are the following (2N + 1) x (2N + 1) pentha- diagonal mass and stiffness matrices 
/ 
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We introduce the discrete analogue of -ffo(0, 1), L 2 (0, 1) and H 1 (0, 1) spaces to be 



Ttu 



0. The inner 



n.^F^f^K^er+'si \Wh\\h,i < 00}, i = -1,0,1. 

For the elements of the space "H^ we have to impose the additional requirement F = F N+1 
products defining the discrete spaces H h , i = — 1,0, 1, are given by 

(Fh,Gh)h,i ■— {{MhSfr 1 ) 1 1 Si l Fh,Gh)i 

i ■ V A ' hi A ' h)h,i: lui < 

Yl'kLi 1 Fk/2Gk/2 (the overline symbol denotes complex 



j C 2N + l , 

and the norms are given by | \F h \\f i :— (F h , F h ) hA , for alii = —1,0,1. Here, ( v ) C 2»+i is the inner product in 

the Euclidian space C n , defined by (F/,, Gh)c 2N+1 
conjugation). 

For //j e Uh with coefficients F^ = {Fj / 2 )i<j<2N+i G we introduce the following notations for the 
three possible discrete derivatives on each nodal point (the forward, the backward and the midpoint one): 

F i -iF j+1/2 +3F j+1 

v h ± j .— jh,x\^j-rj — h _ I . / /, . ■■ .. ... . ' 

d/i-Fj + 1/2 



fh,x{Xj+l-) 



and the values of //, at X7+1/4 

3 



(j + l/4)h and ^+3/4 := (j + 3/4)h, < j < N: 



h 



Fj+l/4 ■= fh{x j + i/A) 



F, 



3 



F 3+3/4 '■= fh{%j+Z/i) 



:F 4 



3 

4 F J + l/2 



3+1- 
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With these notations, it is easy to check that the || • \\h,i and || • H^o - norms admit the following 
representations: 

N 

= I E + 4|5,F J+1/2 | 2 + \d-F J+1 f) 



(12) 



3=0 

N 



M,o = 5& E (7|F,| 2 + 32|F, +1/4 | 2 + 12|^ +1/2 | 2 + 32|F J+3/4 | 2 + 7|F :)+1 | 2 ) 

3=0 



Set Vh '■= T~L\ x H° and its dual V' h :='K h 1 x the duality product (•, -)v' h ,v h between V' h and Vh 
being defined as ((F h ,x, Gh,i), (Fh,2, G h>2 ))v' h ,v h F h , 2 )h,o + (G h ,i, G h , 2 )h,Q- 



Problem (10 1 is well posed in H h x %9. The total energy of its solutions is conserved in time: 



(13) £h(V° h ,ul) ■= l(\\v h (t)\\l A + \\iJ h , t (t)\\io)' 

One of the goals of this paper is to study discrete versions of the observability inequality Q of the form 

T 

(14) £h(V° h M) < °h{T) J \\B h V h (t)\\ 2 c2N+1 dt, 



where By L is a (2iV + 1) x (2N + 1) observability matrix operator. The observability inequality (14) makes 
sense for rather general matrices B^, corresponding, for instance, to the observability from any open subset 
contained in the spatial domain (0, 1). But, within this paper, we focus on the particular case of boundary 
observation operators Bh, in the sense that they approximate the normal derivative u x {x,t) of the solution 
of the continuous adjoint problem ^ at x = 1. The main example of such boundary matrix operators Bh 
that will be used throughout this paper is as follows: 



(15) B 



-i (i s j) = (2JV + l,2JV), 
0, otherwise. 



The operator Bh is also the one used for the finite difference semi-discretization in [13] . Let us remark 
that at the discrete level there are different ways to approximate the normal derivative of the continuous 
solution of |3]). Since B^\Jh{t) is a vector and u x (x,t) is a scalar, the way in which Bh\Jh{t) approximates 



u x (l,t) needs to be further explained. Remark that Bh in (15) is almost a null matrix, excepting the 
penultimate component on the last row. The last component of B^V^it), the only non-trivial one, equals 
to Uh > x( x N+i/2i £)• The consistency analysis shows that if / is a sufficiently regular function and fh is its 
quadratic interpolation, then fh,x(%N+i/2) is a first-order approximation of /^(l). 



We are interested in observability inequalities ( 14 ) in a finite, but sufficiently large observability time, 
say T > T* > 0. In this paper we show that, when working on the whole discrete space Vh X Vh, the 
observability constant Ch{T) blows-up as h — > 0, whatever T > is. One of the main contributions of this 
paper is to design appropriate subspaces Sh C Vh on which the observability constant Ch{T) is uniformly 
bounded as h — > 0. 

We will also prove that the discrete version of |5]) below holds uniformly as h — > in the approximate 
space V h for all T > 0: 

T 

(16) c h (T) J \\B h Uh(t)\\ 2 c2N+1 dt < S h (U° h ,Vl). 



Once the observability problem is well understood, we are in conditions to address the discrete control 
problem. For a discrete control function Vh, we consider the following non-homogeneous problem: 

(17) M h Y hM {t) + S h Y h (t) = -B*V h (t), Y h (Q) = Y°, Y M (0) - Y*. 



Here, the superscript * denotes matrix transposition. Multiplying system (17) by any solution XJh(t) of 



the adjoint problem (10), integrating in time and imposing that at t = T the solution is at rest, i.e. 



(18) ((Y M (T),-Y fe (T)),(U h! U^)) V ;,v, =0, V(U°,U*) e V h , 
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we obtain that ~Vh{t) necessarily satisfies the following identity which fully characterizes all possible exact 
controls Vh(t): 

T 

(19) J (V h (t), J B h U h (*))c^+^t=<(Y^-YO) ! (U ?l (0),U M (0))) v ,, Vh) V(U°,U*) G V h . 
o 

In view of this, we introduce the following discrete version of the quadratic functional ([8]): 

T 

(20) MKM) = \J l|flfcU h (t)||^ +I dt-<CYi,-Y2),(U h (0),U M (0))) vj - iVh) 



XJh(t) being the solution of the adjoint problem ( 10 1 with initial data (U?, Ui) and (Y^, Y?) e the initial 
data to be controlled in (17 1. 

The functional Jh is continuous and strictly convex. Thus, provided it is coercive (which is actually 
what the uniform observability inequality guarantees), it has an unique minimizer (U9,Ui) G Sh whose 
Euler-Lagrange equations are as follows: 

T 

(21) J (B h \J h (t),B h U h (t)) c ™ +1 dt = ((Y{, -Y° h ), (U h (0), U h)t (0)))v; lVjk , 



for all (U°,U^) e Sh and U/,(i) the corresponding solution of (10). The discrete HUM control is then 

(22) V h (t)=V h (t):=B h XJ h (t). 

Let us briefly comment the analogies between the identities ([7| and (22). As we said, when Bh is a 
boundary observability matrix operator, like for example the one in (15), BhXJh(t) is a vector whose last 
component Vfc(i), the only non-trivial one, approximates the normal derivative of the solution to the adjoint 
continuous wave equation |3]). Accordingly, the controls —B^Vh(t) only act on ypf(t) when Vfc(i) is the 
numerical control obtained by ( |22[ ). Consequently, the boundary observability operator Bh does not act 
really at x = 1, but at a; = x./v+i/2j being in fact an internal control acting on a single point which is closer 
and closer to the boundary as the mesh size parameter becomes smaller, so that in the limit as h — > it 
becomes a boundary control. 

Observe also that the control problem we deal with is a coupled system of non-homogeneous ODEs 
modeling the interaction between the nodal and the midpoint components. Thus, the node x N+ i/ 2 lies in 
fact on the boundary of the midpoint component. Consequently, the controls — B^V \{t) in fl7| ), with ~Vh(t) 
as in ( 22 ) , are really natural approximations of the continuous boundary controls v in ([!]) . 

3. Fourier analysis of the P 2 -finite element method 

For the sake of completeness, we recall the spectral analysis of this quadratic finite element method, 
following [10]. The spectral problem associated to the adjoint system ^ is as follows: 

(23) Find {A h ,<p h ) eRxU h such that {(p h , <j>h)Hl = ^hifih, <j>h)^, V<t>h & U h . 

Due to the symmetry and the coercivity of the bi-linear forms generated by the scalar products (•, -) L 2 
and {-,-)h^ we have Ah > 0. Let iph = {&j/2)'\<j<2N+i De ^ ne components of the eigenfunction (ph- The 
pair (Ah,<fh) is a generalized eigensolution corresponding to the pair of matrices (Sh,Mh), i.e. 

(24) S h tp h = A h M h (p h . 



Consider the normalized eigenvalues A := h Ah- The system (24) is a pair of two equations: 

8 - 16- 8. . , 1 . 8 . 1 _ v . 

( 25 ) - 3</>j + y Vi+i/2 " 3^+i - A (^<ft + 15^+1/2 + 15^ +1 ) = °' 

for < j < N, and 

1 . 8 _ 14 _ 8 . 1 _ . / 1 . 1 _ 4^1^ 1 - s 

(26) - Vj _ 1 --y j _ 1/a +— (pj -- ( p j+ y 3 +- Vj+1 -A(-— Vj _ 1 +—^ 
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depending on 1 < j < N and with the boundary condition tpo — <pn+i = 0. 

The acoustic and optic modes. From (25), we obtain that for A 7^ 10, the values of (ph at the 
midpoints can be obtained according to the two neighboring nodal values as follows: 

40 + A 



(27) 



8(10- A) 



(28) 



Replacing ( 27 1 into ( 26 1 , we obtain that 
3A 2 - 104A + 240 



1 . 

2^'" A 2 + 16A + 240 



-(fj+i =0, VI < j < N, with <p = ip N+ i = 0. 



It is easy to check that fij — sin(kTTXj), with 1 < k < N, solves (28 1. Then the normalized eigenvalues 
A verify the identity 

„ » . >, . , 3A 2 -104A + 240 

(29) cos(kirh) = w(A), with w(A) = 



A 2 + 16A + 240 

For each n € [0, n], consider the following second-order algebraic equation in A 
(30) (3 - cos(r/))A 2 - 2A(52 + 8 cos(?/)) + 240(1 - cos(?/)) = 0, 

whose solutions are 



Afa): 



(31) 



A = A a (ri) 



120 sin 2 (77/2) 



and A = A (ry) 



22 + 8 cos 2 (77/2) + 2 V /A(r ? ) 



11 +4 cos 2 (?y/2) + y/A(r]) " ''' l + sin 2 (r7/2) 

with A(ry) := 1 + 268 cos 2 (?7/2) — 44cos 4 (?y/2). The superscripts a and o stand for acoustic and optic. We 
will also need the square roots of the Fourier symbols (31 1, the so called dispersion relations: 



(32) X a ( V ) := \/ A a (n) and X°(r,) := \/A°(rj). 
For each 1 < k < N, set 

(33) A a ' k := A a (knh) and A°' k := A°{k-Kh). 

We refer to (A"' )i<fe<jv and (A°' k )i<k<N as the acoustic and the optic branch of the spectrum. The 
corresponding eigenvectors are 



(34) 



~a.k 



~a,k _ 40+A° 



sin(fc7re,), < j < N + 1, and 'P j ' +1/2 = 4(10 _ A a. )i) 
sin(fcTr^), < j < N + 1, and p>°' k - 40+A °' fc 



cos sin(fc7r.T J+1/2 ), < 3 < N, 



9 



4(10-A°. fc ) C0S (2 



) sin(fc7 



c j + l/2), 



< j < N. 



The resonant mode. Up to this moment, we have explicitly calculated 2N solutions of the eigenvalue 



problem (24 1. To do this, we have supposed that A 7^ 10. But A = A r 



10 is also an eigenvalue 



corresponding to the resonant mode, the superscript r standing for resonant. The components of the 
corresponding eigenvector ip r h are 

(35) <% = 0, < j < N + 1, and <p r j+1/2 = (-1)*, < j < N. 

For any normalized eigenvalue A, define A/, := A/h 2 , X := %/A, := ^/A^ (see Figure [2]). 
Normalized eigenvectors. For any eigenvector (ph g {'pt' > ^>°' k 1 < k < N, <f> r h }, we define the 



L -normalized eigenvector ip^ :— <Ph/\\<Ph\\h,o- 

Using the expression (12 1 of the discrete norm 

\h,o and 



and 



\h,l 



and the identity (27), we obtain the 



following representations of the || • \\h,o and || • H/^i-norms of the acoustic and optic eigenvectors in terms of 
their nodal components, for all a £ {a, 0} and 1 < k < N: 



(36) 
and 
(37) 



I ~a,k 1 12 
\<Ph Wh,0 



1 

24 



1 



500 



110- A' 



,fc|2 



N 
3=0 



~a,k |2 



h 

12 



N 



El ~a,k ~a 



~a,k\2 



3=0 



N ~a,k 



3=0 



^3+1 ^3 
h 



~ot,k „ 

<Pi 2 



1 4 



5A a > k . 

4(10- A a < k ) 



N 

3=0 



~a,k\2 
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Figure 2. The eigenvalues (left) versus their square roots Ah (right): the continuous 
ones (blue), the acoustic (red), the optic (black) and the resonant (green) modes. 



Normalized acoustic eigenvector, k=10 Normalized acoustic eigenvector, k=50 Normalized acoustic eigenvector, k=97 




Figure 3. Examples of normalized acoustic, optic and resonant eigenvectors for N = 99. 
In blue/red, we represent their nodal/midpoint components. 
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Then, using the representation formula ( 36 1 , the identities 

N N 

h sin 2 (kTTXj) = \, h I sm(kicxj + i) - sin(kirxj)\ 2 = 2sin 2 (^p), 
(38) j - Q , v " 

h £ |sin(/c7ra; i+ i) +sin(fc 7 ra; j )| 2 = 2cos 2 (^),V1 < k < N, 

3=0 



and (29 1, we obtain 



-«,*ii2 _ * (A-10)(A 2 + 16A + 240) 



(39) \\0i' \\ho = -^ , withW(A) = 

17 UYh Uh '° 3W(A a ' k ) 19A 2 + 120A - 3600 

Let us remark that ||<^' fc |U,o blows-up as kh — >■ 1. With the above notation, 

(40) ip"' k — n a,k sm(kirxj) and ^"+1/2 = mQ,fc sin(fc7rXj +1 / 2 ), Va € {a,o}, 
where 

n and m standing for the nodal and midpoint components. 

Using the explicit form of the || • \ - norm and the characterization of the resonant mode (35), we 
obtain that ||<^|| 2 = 8/15 and therefore, the normalized resonant mode ip r satisfies 

(41) ip r 3 = 0, V0 < j < N + 1 and <p r j+1/2 = (-l) j \/l5/2\/2, V0 < j < N. 
Let us introduce the sets of eigenvalues and of 7^-normalized eigenf unctions, i.e. 

(42) £V h := {A^ k , a € {a, o}, 1 < k < N, A r h } and £ T h := K' fe , a e {a, o}, 1 < k < N, cp r h }. 

Fourier representation of discrete solutions. Since SJ-h is an orthonormal basis in "H", the initial 
data in ( |10[ ) admit the following Fourier representation: 

V N 

(43) V\ = u a ' k > i <p a h' k + V'Vh + Y, 2°'*'Vfc' fc , where = (U* h , ^' fc ) M and vT>* = (U* h , ^) M , 

fc=i fe=i 

with a € {a, o}, 1 < fe < N and i = 0, 1. Correspondingly, the solution of ( |To| can be represented as follows: 

V N 

(44) U h (t) = ^ [^M± fc exp(±ia°'' £ )^ fe +i4exp(±itA^)^ + ^ exp^A^ 

± fc=i fc=i 

where 

1 



The total energy ( j 13[ ) of the solutions of ( 10 ) is then as follows: 

V N 

(45) £ h (u° h M) = E A £'Wl 2 + ) + K(\K\ 2 + \u-\ 2 ) + E A f (K/l 2 + l^ fe l 2 )- 
fe=i fe=i 

Limits and monotonicity of the eigenvalues. Firstly, let us remark that as kh — > 1, A a,fc — » 10, 
A°' k — > 12 and as fc/i — > 0, A° ,k — > 60. On the other hand, the so-called group velocities, which are first-order 
derivatives of the Fourier symbols (31 ) or of the corresponding dispersion relations (32 ), verify the following 
positivity condition 

d v \ a (v), d v A a ( v ), ~d v \°( V ), -^A°(77) > 0, V77 e (0,tt), 
which means that the acoustic branch is strictly increasing and the optic one is strictly decreasing in k. Con- 
sequently, the high frequency wave packets involving only the acoustic or the optic modes and concentrated 
around a given frequency 1 < k* < N propagate in opposite directions. Moreover, at rj = or rj = it, the 
group velocities satisfy 

9„A a (7r) = d v A°(7r) = <9 r ,A a (0) = <9„A°(0) = 0, <9„A a (7r) = d v \°(n) = 9^A°(0) = and d v X a (0) = 1, 
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which, according to the analysis in |17) . shows, in particular, that there are waves concentrated on each 
mode which propagate at arbitrarily slow velocity. 

4. Boundary observability of eigenvectors 
The main result of this section is as follows: 



Proposition 4.1. For all a 6 {a,o} and all 1 < k < N, the following identity holds for both acoustic and 
optic eigensolutions: 

I <Pn 2 



(46) ||^' fc ||L = 



W{A a > k ) 



with W(A) = 



24(A-10) 2 (A-12)(A-60) 
(-19A 2 - 120A + 3600)(A 2 + 16A + 240) ' 



Moreover, for the resonant mode, the following identity holds: 



(47) 



lG 



Vn+1/2 



Remark 1. The identity (46) is the discrete analogue of the continuous one \\<f k \\ 2 H i — \(p k (l)\ 2 /2, where 
ip k (x) = \/2sm(kTTx) is the L 2 -normalized eigenfunction corresponding to the eigenvalue A k = k 2 ir 2 . 

Remark 2. Due to the monotonicity of the Fourier symbols, we have that A a ' k £ (0, 10) and A° ,k £ (12, 60), 
for alll <k< N. The quadratic equation -I9x 2 - 120a; + 3600 = has the roots x x = -60(1 + \/20)/19 < 
and x 2 = 60/(1 + 720) € (10, 12). This allows us to guarantee that W(A) > 0, for all A £ (0, 10) U (12,60). 



Remark 3. Due to the form of the denominator in the right hand side of \4ty and from the above lower and 
upper bound of the Fourier symbols and the behavior of the group velocities, we deduce that the coefficient of 
\ip°^ k /h\ 2 in the right hand side of (4.6) is singular as kh 1 both when a = a or a = o and when kh — > 
and a = o. 



Proof of Proposition 4-1 Fix a £ {a, o} and 1 < k < N. Obviously, it is enoug h to prove (46) for the un- 
normalized eigenvectors ip a,k . We will use two approaches to prove the identity (46). The first one consists 
on using the classical multiplier Xj — (pj-i) /2h (which is a discrete version of the continuous one Xip x ) 
in the simplified spectral problem (|30|) and then to apply the Abel summation by parts formula 



(48) 



N N 

y^aj+i - a,j)bj = a N+1 b N+1 - ai&o - aj+i (fy+i - bj), 

3=1 3=0 



for all (aj)i<j<jv+i £ C N+1 and (bj)o<j<N+i £ C N+2 . In what follows, we will add the superscript a, k to 
the solution ip^ of ( 28 ) . In this way, we deduce the following identity: 



iV 



3=0 



1 I ~a,k\ 2 
2 1 " 



iV 



<py\' + (w{K* k )±i)hY i P i &<pT k - 



3=0 



Replacing the crossed sum h^2^ =0 <pj£ 1 ipj' k obtained from the identity with + into the one with 
the following equality: 

(49) 



we get 



3=0 



N ~a.k 



<Pi 2 



w(A a ' k ) + 1 



1 w(A a ' k ) 



1 N 

1, \ ^ I ~a.k 

h?w(A^ k ) + i n 2^m+i 

^ ' 3=0 



~a,k I 2 



Using the representation (37) of the || • H^^-norm of the optic and acoustic eigenvectors in terms of the 



nodal components, we obtain that (|49|) is equivalent to 
(50) \\<P% k \\Li 



w(A a > k ) + 1 



~OL,k 

<£n_ 
h 



1 /w(A"' fc ) - 1 
h?\w{A a ' k ) + 1 



5A a,fc 

4(10- A a - k ) 



N 

7 I ~ot,k 

hy \ip,' 



3=0 



'3 + 1 



~a,k I 2 

■<Pi • 
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Replacing the representation of the 1 1 • | j/i^-norm of the eigenvectors ( 36 ) into the one of the 1 1 • | j/^i-norm, 



(37), we obtain: 



II -a,fei|2 l2 u~^ k \\ 2 1 25A"- fc (A^ fc - 12) - 6(A a ' k - 10)(A^ fc - 60) 

(5i) IK || M = j^Wr \\ h ,o + p 12(AQ , fc - 10)2 



V 



L \ ~* I ~a,k 



~a,k I 2 

■<Pi ■ 



3=0 



On the other hand, the || • H^o and || • ||h i-norms of the eigenvetors are related as follows: 

A a > k , 



(52) 



I -ot.k 1 1 2 
\<Ph ll/i.l 



h 2 



I ~a,k 1 12 
\fh Wh.O- 



Replacing (52 1 into (51 1, we get 



N 



(53) h^T \(pf^ +tf' k \ 2 



3=0 



12(A Q < fc - 10) 2 (A"' fc - 12) 
A a > k [-6(A a - k - 60)(A Q < fe - 10) + 25A Q < fc (A Q ' fe - 12)] 



t 2 1 1 ~a,k 1 12 

h \\<p h ' 



By combining (|50|) and (531, we obtain 



1 



6(A Q > fc - 60)(A a > fc - 10) + 25A Q ' fc (A. Q ' fc - 12) 



-6(A Q . fc - 60)(A Q < fc - 10) + 25A«' fe (A Q ^ - 12) 
from which the identity Mm follows immediately. 



I ~a,k 1 12 
\<Ph \\h,l 



(A a - fe ) 2 + 16A Q > fc + 240 



2(A«. fc - 10)(A«< fe - 12) 



~a,k 



the relation (29 1. Thus, for W as in (|46|), we get 



The second approach to prove (46) is much more direct. It c ons ists in using the representation (371 of 
the || • ||/ l; i-norm of the eigenvectors, the trigonometric identities (38), the fact that \{p°^ \ = | sin(fc7i7j)| and 



-«,fc||2 l-w(A"' fc ) + | 



5A°' fc 
4(10-A°. fc ) 



'(1 + w(A"' fe )) 



I 'Pn 



(1 - w(A a > k ))(l +w(A a > k )) 



W{A a ' k )' 



The identity ( 47 ) follows by combining the explicit expressions of the components of the resonant eigen- 

□ 



vector (35) and (12). This concludes the proof of (47) 



5. Discrete observability inequality: an Ingham approach 



The observability inequality. In this section we prove that the discrete observability inequality ( 14 ) 



holds uniformly as h — » in a truncated class of initial data for the observation operator Bh introduced in 



(15). More precisely, consider < A3. < 10 and 12 < A° < A+ < 60 and correspondingly the wave numbers 



(54) n% := {A a )-\A a + ) and V ° ± := (A )" 1 ^) 
and introduce the subspace of <C 2N+1 given by 

%,n a ,vl,vl ■=spa,n{ip a ' k ,kirh < ?7 + } © span{( / 5°' fc , 77° < kirh < rf_ }. 
Consider the truncated subspace Sh C Vh defined by (see Figure [4j: 

(55) Sh ■= (Th^,ri _,-q + X Th,7)^,r] _,7)^) H Vh- 

The main result of this subsection is as follows: 

Theorem 5.1. For all A\ e (0,10) and A° < A% G (12,60) independent of h, all initial data (\5 h '° ,\J h ^) e 



Sh introduced in (55), with 77° , 77^ given by (54), and all 



n+,V-,V- mm { mm d v X a , min {-d v X°(r]))} 



the observability inequality \ 14) holds uniformly as h — > for the operator Bh in 1 15) 
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Figure 4. The selected area contains square roots of eigenvalues whose eigenvectors gen- 
erate the truncated class lh,rf ,nl-v° ■ ^ n red/black/green, the acoustic/optic/resonant mode. 



Proof of Theorem \5.1\ The fact that the eigenmodes involved in the class Th^,ri o ,77° & r e such that the gap 
in each branch has a strictly positive lower bound uniformly as h — > allows us to apply Ingham Theorem 
(cf. |15j . Theorem 4.3, pp. 59). More precisely, the spectral gap on each branch is bounded as follows: 

l "' H1 -Af >tt min d n X a {n) > 0, Vkwh, (k + l)nh € [0,r£] 



A',' 



ne[o,vi 



and 



\0,k \ 

A h - A h 



o,fc+l 



> 7T min {-d n X°(rj)) > 0, Vfor/i, (k + l)nh G [rf + , rf_ 

J7S[t)°,»;°] 



Then 7 = 7(77°, 77° , 77° ) := 7rmin{ min d v X a {rj), min (— d v X (r]))} > is the uniform gap needed 

»7e[o,j?+] v€.[n+, , nl] 

to apply the Ingham theory. From (|44| and the definition of the class Th.r^ ,77° ,77° > we have 



E 



~-o.k 



/ 1 - . \ o,k \ o,k 



=t k-KhG [0,77" ] kirhG [77^,77° 



By applying the inverse inequality in Ingham Theorem (cf. [TS], pp. 60, (4.9)), we can guarantee that for 
all T > 2ir/j — T* a v o ^ , there exists a constant C_(T) > independent of h such that 



tf-cn£[ E 



^a,fc i2 



a.k 

<Pn 



h 



E > ?i 

fc7r/tS [77^,77° ] 



^o,fc i2 



o,k 



h 



< 



U N (t) 



dt. 



Using the identities (46) and H^' ||^ x — AJ*' , we obtain that (the function W below is as in (46 )): 

T 



± fc7r/i<E[0,7;° ] fcT^e [77^,771] 



°' fc l^ fe | 2 ^(A°^) 



<-I 


U N (t) 




h 



dt. 



Taking into account that for our filtering algorithm A E [0, A+]U[A° , A+] c [0, 10)U(12, 60) (with uniform 
inclusion as h — > due to the fact that A^ and Ajj. do not depend on h), and that for A G [0, 10) U (12, 60), 
the weight W is strictly positive, we can guarantee that 



AS[0,A» )U(A° ,A° ) 



W(A) is strictly positive. 



Then the proof concludes by taking in (14) with Bh given by (15) the observability constant C/j(T) 
C_(T)C_ (77+, 77° , 77° ) which is independent of h. 



□ 



The admissibility inequality. Using the direct inequality in Ingham Theorem (cf. 15J, pp. 60, (4.8)), 
we can also prove that the inequality (16), with Bh as in (15), holds uniformly as h — > for all T > 0, with 
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l/c/ l (T) = C+(T)C+ (77° , 77° , 77° ), where C+(T) is the constant of the direct -Ingham inequality and 

C + (r)l,ril,r)i) := max W(A) is bounded. 

+ + Ae[0,A^)U(A° ,A°) 

6. Discrete observability inequality: a bi-grid algorithm 

1. The observability inequality. In this section, N will be an odd number. We consider the space 
Ch containing piecewise linear functions given below 

C h := {F h = 0F i/2 )i<i<2W+ 1 with F = F N+1 = 0, s.t. F j+1/2 = (Fj + F j+1 )/2, V0 < j < N} 

and the space Bh of the discrete functions whose nodal components are given by a bi-grid algorithm, i.e. the 
even components are arbitrarily given and the odd ones are computed as average of the two even neighboring 
values: 

B h := {F h = (F J/2 ) 1 < J < 2JV+1 with F = F N+1 = 0, s.t. F 2j+1 = {F 2j + F 2j+2 )/2, V0 < j < (N - l)/2}. 
We also define the subspace Sh of Vh 

(56) S h := ((C h n B h ) x (C h n B h )) n v h . 



The aim of this section is to prove that the observability inequality ( 14 1 still holds uniformly as h —> 



for initial data in the bi-grid subspace Sh introduced in ( 56 1 : 

Theorem 6.1. For all T > 2 and all initial data (U^Ujj) in the adjoint problem (10) belonging to Sh 



introduced in (56), the observability inequality (lj.) with Bh given by (15) holds uniformly as h — > 0. 

Remark 4. Note that, in the bi-grid filtering mechanism we have designed, the data under consideration 
have been taken, before filtering through the classical bi-grid algorithm, to be piecewise linear in each interval 
(xj,Xj+\), j € Z, which imposes a further restriction. This allows to obtain the sharp observability time. 

The bi-grid filtering algorithm proposed in Theorem \6.1\ yields optimal observability time, i.e. the char- 
acteristic one T* = 2. This is due to the fact that for a numerical scheme the minimal time required for 
the observability to hold is 2/v, where v is the minimal group velocity involved in the corresponding solution. 
From our analysis, we will see that the bi-grid filtering algorithm above acts mainly as a Fourier truncation 
of the whole optic diagram and of the second half (the high frequency one) of the acoustic one. Consequently, 
v := min^ e [o l7r /2] d^X 1 ^). Since the group velocity of the acoustic branch, d v X a (r]), is increasing on [0,7r/2], 
we conclude that v = d v X a (0) — 1 and then the observability time of the numerical scheme is sharp: T > 2. 

The following two auxiliary results hold: 



Proposition 6.2. If the initial data (U?,Ul) in (10) belong to {Ch x Ch) HVi,, then the resonant Fourier 



coefficients in (44) vanish, i.e. 

(57) tf ± = 

and the optic and acoustic ones are related by the following two identities: 

(58) (u^ + 21 ,fc )(m a ' fc - 7^< fc cos (^)) + (uf + v°_! k )(m°> k - n°- k cos (^)) = 



(/// 



d 



(59) \ a < k (v% k - u a _: k )(m a ' k - n a ' k cos (*^)) + X°' k (u°( k - v°_l k )(m°> k - n°< k cos (^)) 



(60) 



Taking squares in (58 1 and (59) and in view of (401, we deduce that 

-a,fci2 __j i-o,fc -o,fe,2 _ A a,fc Wi(A a ' k ) ! 



\u, + u_ 



where 



Wt(A a > k ) k 

vKi(A°. fe r + 



and \u + 



-~a,k\ 2 



A°' k VKi(A°> fc ) 



A 2 fA 2 + 16A + 240) 
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Proof of Proposition 6.2 We will prove only (57) and (58), the proof of (59) being similar to the one for 
(58). Observe that the Fourier representation of the identity that characterizes U° £ Ch is 

N 

(61) = ^[(^+ fe + u a _i k ){m a ' k - n a ' k cos{knh/2)) + (v% k + u°! k ){m°> k - n°' k cos(/ctt/i/2))] sin(fc7ra; i+1/2 ) 
\/l5 

^ 271 sin (( N + 1 ) 7rx J+i/2)' 



k=l 



for all < j < N. Multiplying (61) by sm(lTTX j+1/2 ), 1 < I < N + 1, adding in < j < N and taking 



into account that h^2^_ sm (^ 7ra; j+i/2) sin(Z7riEj +1 / 2 ) — ^fc,i/2, for all 1 < k, I < N + 1, we conclude the two 
identities d57| and {f58|) . □ 



The total energy of data (U2,Ujj) € x £/j) n Vh in (10) can be written only in terms of the nodal 
components and coincides with the one of the Pi-finite element method 



(62) 



u N ,TfO tt° 



h 
12 



N 



J2(W}\ 2 + \u} +1 + u}\ 2 ). 



3=0 



Taking into account the form of the Fourier coefficients (57) and (60) corresponding to linear initial data, 



we obtain that the Fourier representation of the total energy (62) is as follows 

N 



(63) 



* t DS,u!)-i£Ar'[(i + £K)W 



k . ^-a,k\2 
+ U' 



1 



Wi(A a -' 



i^a,k ^a.k\2 



The second auxiliary result establishes that for initial data in the bi-grid subspace Sh in (56), the high 
frequency Fourier coefficients on the acoustic branch can be evaluated in terms of the low frequency ones: 



Proposition 6.3. For each element (U°,U^) € Sh introduced in (56), the following identities hold: 
(64) 



r (N+l)/2 = r JN+l)/2 = 



(65) 
and 
(66) 



^a,N+l— k , ^a,N+l-k 



^.N+l-k VK 2 (A a ' W+1 - fc ) k 



(u+ + uj ), 



i a ' k W 2 (A a ' fc ) 



„a,JV+l-fc u/ I \a,N+l-k\ 



for alll<k<(N- l)/2, with 



W 2 (A) 



n a > k W 2 (A a ' k ) 

(60-A)(A-10)(A-12) 
(A 2 + 16A + 240) 2 



Taking squares in ( [651661 ), we obtain that 



(67) 
and 



f ao \ \^a,N+l-k ^~a,N+l-k\2 

(68) |m + — Uj I 



\^a,N+l-k , ~a,N+l-k\2 W 3 (A a - N+1 k ) a k ~a,k\2 



1 W 3 (A a > k ) 

A a,k W3 ( A a,N+l- k) 



^a,k\2 



I 2 , VI < k < (N — l)/2, 



where 



W 3 (A) = Wi(A) 



I fct,N+l-k W 3 (A a ' k ) ' + 

(A - 10)(A 2 + 16A + 240) _ (60 - A) 2 (A - 10) 3 (A - 12) 2 



19A 2 + 120A - 3600 (19A 2 + 120A - 3600) (A 2 + 16A + 240) 3 
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Proof of Proposition 6.3 Due to the orthogonality in °f the eigenvectors in £Th, to the fact that G Ch 
and using the representation (40) of the normalized eigenvectors, the following identity holds 



N 



(69) 



tO a,k' 



^a,k , ^a,k / TT 0,fc\ 4(60 A a ' ) a/ciV^rrO • n \ 

l + + u - = ( U W )h,o= ( A a, fc) 2 + 16A a, fc + 26 o n h^U^smiknx,). 

3=1 



Now, using the fact that U° € Bh, we have 



v 



h Uj sin(knxj) = 2 cos 
3=1 



(A-l)/2 

J ' f) h E ^W^o 

3=1 



Taking (29) into account, we have 
(70) 



*^a,k . ^a,k 
U, + U_ 



(N-l)/2 

l6W 2 (A a ' k )n a ' k h J2 U^sin(knx 2j ). 

3=1 



For k = (N + l)/2, we obtain u a _; (N+1)/2 + u"' (JV+1)/2 = 0. Since sin((7V + 1 - k)TTx 2j ) = -sm(kwx 2 j), we 
obtain (65) by equating the expressions of ^2i'=i U®j sm(kTTX2j) form the identity (70) corresponding 
to k and to N + 1 — k. The proof of (66 1 is similar, based on the fact that (U^, (p a,k )h,o — i^' k {u+ k — u°2 k ), 



from which for k — (N + l)/2 we obtain that u^ N+1 ^ 2 — ^( A '+ 1 )/ 2 _ w hj cn concludes (64). 



□ 



Replacing the Fourier representations ( 67 ) and ( 68 ) into the total energy with linear initial data ( 63 1 , 
we obtain that energy corresponding to initial data (U^U^) € Sh in (56) is given by 

(W-l)/2 

fc=l 



where the low and high frequency coefficients are given by 

) wlo = W^) 
)' rv -> k - + VFi(A°^)' 

^a,N+l-k ^ fro,N+l-k Wi(A a > N+1 - k )\ W 3 (A a ' N+1 - k ) 



,,. , , A^W^A^) w;o __ i , W l {A a ' k ) wM _^ , W 1 (A a - N+1 - k ) \ W 3 (A a ' N+1 - k ) 

W 1 (A°> N + 1 - k ) J W 3 {A a ' k ) 



+> k A a k W^A"^) 



W1\ := 1 



and 



Jya,k \ j±a,N+l-kwJAo,N+l-k\) W 3 (A a > k ) 



For any 5 € (0, 1) which does not depend on h and any solution (44 ) of ( 10 ), let us introduce its projection 
on the first 5(N + 1) frequencies of the acoustic branch to be 



(72) 



<5(iV+l) 

r a s u h (t) = J2 E u a /cM±*t\ a h k )v a h k 

± k=l 



Observe that the projection (72) is still a solution of (10), therefore its total energy is conserved in time. Set 



<5(JV+1) 



4(r^(U°,Ui)):=f h (rfU ft (0),r^U M (0))= £ A a h k (\u a + 



k=l 



The following result provides a bound of the total energy of the solutions of ( 10 ) with initial data (U^, U^) € 
Sh as in (56) in terms of the total energy of their projections on the first half of the acoustic mode: 



Proposition 6.4. For any solution XJ h (t) of (10) with initial data (U",U^) E Sh introduced in (56), there 
exists a constant C > which does not depend on h such that 



(73) 



^(u^u^c^r^u",^)). 
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Proof of Proposition 6.4 In order to obtain the estimate (73), we claim that it is sufficient to prove that 
there exist W£, > independent of h such that W£ k < and Wf M < for all 1 < k < (N- 1)/2. 
Assuming this for a moment, we can take C := maxjW 7 ]?, W±} for which (731 holds. 

Let us analyze the boundedness of the four coefficients. Observe that W_ k involves the product of 
Wi(A°' fe ) and l/W x {A°' k ) for kh < 1/2. But Wi(A) is singular only for A -> 10, whereas for kh < 1/2, due 
to the increasing nature of A a , we have A"** < A a (7r/2) = 60/(13 + 2^31) < 3. Also l/Wi(A) is singular as 
A -> 0, but since A° is decreasing in 77 we have A°- k > A (vr/2) = (52 + 8\/31)/3 > 30 for all kh < 1/2. In 
the same way, W^ k contains the product of W 1 {A a - k ) / A a - k with A°< k /Wi(A°< k ), for kh < 1/2. The second 
factor A/Wi(A) has a singularity as A — ► 0, but is evaluated for A = A° ,k which is far from the singularity 
for all kh < 1/2. Since VFi(A) contains a factor A 2 , we deduce that Wi(A)/A has the same singularities 
as Wi(A), i.e. A = 10, but since we work on the first half of the acoustic diagram, A Q:fe is far from that 
singularity. We conclude the existence of W l ± > independent of h such that W± k < W±. 

The coefficient W_ k contains two terms. The first of them is constituted by the factors W 3 (A a ' N+1 ~ k ) 
and l/W 3 {A a ' k ). Since V 3 (A) is not singular for any A G (0,10) U (12,60), then W 3 {A a ' N+1 ~ k ) is bounded 
for all kh < 1/2. On the other hand, l/W 3 (A) has three singularities: A = 10, A = 12 and A = 60. Since 
all A-s situated on the first half of the acoustic mode are far for all the three singularities, we deduce the 
boundedness of l/W 3 (A a ' fc ). The second term is a product of four factors: Wi(A a ' N+1 ~ k ), W 3 (A°- N+1 ~ k ), 
l/Wi(A°- N+1 ~ k ) and l/W 3 (A a ' k ). In view of our previous analysis, we deduce the boundedness of the last 
three factors. The first factor blows-up like 1/(10 — A a ' 7V+1 ~' c ) for small k. Nevertheless, in the same range 
of k-s, W 3 (A a,N+1 ~ k ) behaves like (10 — A Q:Ar+1 ~ fc ) 3 , compensating in this way the singularity of the first 
factor W 1 {A a ' N+1 - k ), so that W 1 {A a ' N+1 - k )W 3 (A a > N+1 - k ) is bounded for all kh < 1/2. 

The coefficient W^ l k also contains two terms. The first of them includes the factors: \a-M+i-k^ 
l/A a ' k , W 3 {A a - N+1 ~ k ) and l/I^ 3 (A Q ' fc ). We have already analyzed the first, the third and the fourth 
ones. The second one blows-up like sm~ 2 (knh/2) for small k. But, as we said, W 3 (A a - N+1 ~ k ) behaves 
like (10 — \a,N+i-k^3 ^ s [ n 6 (^knh/2) for small k, compensating the singularity of the second factor, so that 
w ^ A aM+i-ky A a,k ig bounded for all kh < 1/2. The second term in W\\ k contains the factors 1/A°' fe , 
A° ,Ar+1 ~ fe , W 3 (A a,N+1 ~ k ) and l/W 3 (A a ' k ) and is bounded by the same arguments we used for the first term. 

< W±, which concludes the proof. □ 



Consequently, there exist W± 



hl > such that W^ l k 



Remark 5. 



Set CI 



{F h = ( J F J y 2 )i< i <2iv+i s.t. F j+1/2 = a{Fj + F j+1 ), < j < N} and 5^ := n 



Bh) x n Bh)) n V/j. Observe that Lh = C]/ 2 and Sh = Sj/ 2 , where Sh is the space in (56). We want to 



point out that the result of Theorem 6.1 is not longer true when replace Sh by 5?, with a 7^ 1/2, so that the 



condition on the initial data to be linear is sharp. Indeed, when replacing Sh by Sff, in particular Wv? k 



71 ) has to be substituted by 



W 



+,k 



1 + 



A o,fc W ^ K a,k^ 

A a ' k Wf (A°>) ! 



with W?(A) 



1 (A 2 + 16A + 240) (40 - 80a + (1 + 8a)A) 2 



25 



(A - 10)(19A 2 + 120A - 3600) 



One can show that, for a ^ 1/2, it is not longer true that Wi(A a ' k ) — > as kh — > 0, so that this factor 
cannot compensate the singularity of 1/A a,k as kh — > 0, like for a = 1/2. Consequently, for a ^ 1/2, at 
least W?' k is not bounded for 1 < k < (N — l)/2. 



Proof of Theorem 6.1 Proposition 6.3 ensures that the total energy of initial data in Sh introduced in (56) 



is uniformly bounded by the energy of their projection on the first half of the acoustic mode. On the other 



hand, Theorem 5.1 guarantees that the observability inequality (14) with Bh as in (15) holds uniformly as 
h — > in the class of truncated data lying on the first half of the acoustic mode. Combining these two facts, 
one can apply a dyadic decomposition argument as in 12] and conclude the proof of Theorem 6.1 □ 



2. The admissibility inequality. In the rest of this section, our aim is to prove the direct inequality 
(16) for the boundary operator Bh in (15) using the spectral identities (46) and (47). Firstly, let us observe 
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that for all matrix operator Bh and all T > 0, we have the following inequality: 

T T 



(74) 



B h V h (t)\\ 2 c2N+1 dt = / ((B h S^ 1/2 yB h S- 1/2 Sl /2 V h (t),Sl /2 V h (t)) c2 N + , dt 



< \\(B h S- 1/2 )*B h S- 1/2 \y» + t [(S h V h (t),lJ h (t)) c2N+ tdt 



< 2T£ h (Vl,V 1 h )\\{B h S- 1/2 yB h S h 



-1/2, 



|C2JV + 1. 



For any matrix A, the matrix norm || • ||c2n+i involved in the right hand side of (74) is defined as 

1/2 

||A* J 4|| C 2w+i = max HAM,/ iph\\c 2N + 1 - 

ll¥>hlU,o = l 

In the above definition of the norm of a matrix, we can reduce the set of test functions to iph £ £ Th introduced 
in (42), which is an orthonormal basis in M 2Ar+1 . Let us remark that for any eigenfunction iph £ £J~h, the 

]/2 1/2 1/2 

corresponding eigenvalue A/j £ £Vh verifies the identity S h M h iph = A h iph- Consequently, for any 
matrix Bh and (fh £ £J-h, we have 

(75) \\B h S^ 1/2 MjJ 2 Lp h \\ C 2N+i = K^^WBh^hWciN+i. 



For Bh in (15), using the identities (75), (46) and (47), we conclude that 

2 



WiBhS-^TBhS^Tc^ 



1/2.12 



1 

h^SFh Ah 



{ max W(A),—\ 
L Ae(o,io)u(i2,60) 16 J 



is a quantity independent of h. 



7. Convergence of the discrete controls 



In this section, we describe the algorithm of constructing the discrete controls of minimal L 2 (0, T)-norm 
and we prove their convergence towards the continuous HUM boundary control v(t) in ^ as h — > 0, under 
the hypothesis that both inverse and direct inequalities ( 14 ) and ( 16 1 hold uniformly as h — > 0. As we saw in 



the previous sections, the above hypothesis holds when the initial data in the adjoint system ( 10 ) is filtered 
through a Fourier truncation or a bi-grid algorithm. 



7.1. Description of the algorithm. Using the admissibility inequality (16) and the observability one (14), 



one can prove the continuity and the uniform coercivity of Jh defined by (20) on Sh C Vh, where Sh can 
be both the truncated space (55) or the bi-grid one (56). Moreover, it is strictly convex. Applying the 
direct method for the calculus of variations (DMCV) (cf. [3]), one can guarantee the existence of an unique 
minimizcr (U^,U^) of Jh, i.e.: 



(76) 



'-h ■- 



Moreover, the Euler-Lagrange equation (21) associated to Jh characterizes the optimal control (|22) 

JTo 



Remark that when the space of initial data in |10j) is restricted to a subspace Sh C Vh, for example the 
ones given by (55) or (56), the exact controllability condition (18) holds for all (U°,U^) £ Sh- This does 
not imply that the final state (Y^ t (T), — Y^(T)) in the controlled problem (17) is exactly controllable to 
the rest, but its orthogonal projection T$ h from Vh on the subspace Sh, i.e. 



r Sh (Y h|1 (r),-Y / ,(T)) - 0, meaning that ((Y hit (T),-Y h (T)), (U° h ,Vl)) K y h = 0, V(U£,U£) £ S h . 
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7.2. The convergence result. Set Wft(t) to be the last component of the control V/j(t) in (22 1 (the other 
ones vanish). Since Ih < J^(0,0) = and taking into account the inverse inequality |l4"|), we obtain: 



(77) 



1 

J \v h (t)\ 2 dt<8C(T)\\(Yl-Y] 



h)\\v( > 



where C(T) = sup Ch(T) and Ch(T) is the observability constant in (14 1 under filtering. 

he(o,i) 

Set £ J- := {ip k (x) = \/2 sin(/c7rx)}, A k :— fc 2 7r 2 , A fe :— kir and I 2 be the space of square summable 
sequences. Since SF is a Hilbertian basis in each H S (Q, 1), s € K, the initial data (y 1 ,—y°) G V to be 
controlled in the continuous problem admits the following Fourier decomposition: 

oo 

(78) !/*(*) = £v*V(aO, with^ = (y\/) L2 , i = 0, 1. 

fe=i 

Moreover, their || • ||y'-norm has the following Fourier representation: 

J>°_ /\^k,l\2 

(79) 



\\(y\-y°)\& = Yl 



k=l 



\!J" 



T*,0|2 



A fc 



Since the set £ Fh introduced in (42 1 is a basis in IR 2Ar+1 , the initial data (Yi, Y2) g to be controlled in 



the discrete problem (17) admit the following decomposition 

N N 
/ on \ A/ - \ ^a,k,i a,k . ^rA r , \ ^o,k,i o,fc w- n 

(80) Y fc =2^Vh <Ph + Vh fh + 2^ ^ ^ ' Vz = °' X ' 



fc=l 



with 



= (Yi,<p% K ) h , , aeM,l<K iV, and = (Y£, ^) M , Vi = 0,1 
Moreover, their || • ||y<-norm can be written in terms of the Fourier coefficients as follows: 

N 

,,n,9 . 
Wh 



(81) ||(Yi,-Y°)|| v , =J2 lVh 



A 

fe=l 



^r,l|2 I — o,fc, 1 1 2 N 



^o,fc,0|2 



AT 



fe=i 



fe=i 



k=l 



The main result of this section is as follows: 
Theorem 7.1. In the controlled problem [PA), we consider initial data with the following two properties: 



(82) 
and 
(83) 



' ~a,k,l 
\a,k 



k,l 



and (yh' k '°)k — v (y k '°)k as h — > 0, weakly in £ 2 , 



^o.k, 1 



arte? (y^' fc '°)fc as ft — > 0, weakly in I 2 



Then 
(84) 



U/j -^i) as 



0, weakly in L 2 (0,T), 



where v% is the last component of the discrete optimal control Vfc(t) given by (22) obtained by the minimiza- 
tion of the functional Jh on the subspace Sh defined in (55) or in (56) and v(t) is the continuous HUM 
control |[?|). 



Moreover, if the convergences in (82) and (83) are strong, then the convergence of controls in (84) is 
also strong. 
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Proof of Theorem \7.1\ Firstly, let us observe that from (82) and (83), we obtain that there exists a constant 
C > independent of h such that 



(85) 



(Yi,-Y°)|| n <C. 



By combining ( 77 1 and ( 85 1 , we get the uniform boundedness as h 
so that 



of the discrete control Vh in L 2 (0, T), 



(86) v h ^v* as /i->0, weakly in L 2 (0,T). 

It is sufficient to prove that the weak limit v* coincides with the continuous HUM control v given by ([7]). 

The control v can be characterized as the unique control v in ([!]) which can be expressed as the space 
derivative of a solution of the adjoint problem ^ evaluated at x = 1. Then, we have to prove that v* is an 
admissible control of the continuous wave equation, i.e. it verifies the identity ([6]), and that v* = u*(l,t), 
where u*(x,t) is the solution of the adjoint problem Q for some initial data (€,*<", u*' 1 ) £ V. 
Step I - The weak limit v* is a control in the continuous problem ([T]). 

, k £ N} is an orthonormal basis in V := Hg x L 2 (0, 1), then the fact that v is a control 



Since {(±ip /iX 
in 

k £ N. The solution of the adjoint problem (j3j) with this kind of initial data is u{x,t) — ± exp(±iX k (t 
T))ip k (x)/iX k . The condition ^ is equivalent to the following one: 



lj), so it satisfies ^jj, it is equivalent to prove ^ for all initial data of the form (u°, u 1 ) — (±ip k /iX k , ip k ) 



(87) 



v(t)exp(±itX k )dt = 



(-1)* 

V2 \X k 



A:,0 



,Vfc £ N. 



Let us check that v* satisfies (87). To do it, we distinguish between the two cases of subspaces Sh of filtered 



data. Firstly, let us consider the case of truncated initial data, i.e. Sh is given by (55). A particular class of 
initial data in Sh is (U^UjJ = (±y>?' /iA?' , <p^ k ), for which the solution of the discrete adjoint problem 
(10) is U h (t) = ±(p a h ' k exp(±iX a h ' k (t-T))/iX a h ' k , for all kwh < From (21), we see that the discrete control 



Vh verifies the identity 



dh(t) exp(±itA^' fe ) dt = 



(-l) k X a h 



TW h ,Vkirh < rj + 



X 



h 



Let us fix k £ N (independent of h). In that case, exp(±i£A^' fe ) — > exp(±^A fc ) as h — > strongly in L 2 (0,T), 



so that, taking into account the weak convergence (|86j), we can pass to the limit as h — > in the left hand 
side of (88) and we obtain the left hand side of (87) with v substituted by v*. On the other hand, taking into 
account the condition ( 82 ) , which is valid for all test sequences in £ 2 and in particular for the basis functions 
of £ 2 , e k — (0, • • • , 0, 1, 0, • • • ) (meaning that the weak convergence in £ 2 is a pointwise convergence), and 
additionally the fact that n a,k — > y/2 and lXt k / sm(kirh)/h 1 as h for each fixed k, passing to the 
limit as h — >• in the right hand side of (21), we obtain the right hand side of (87). Then, the weak limit 
v* of the optimal control Vh obtained by minimizing the functional Jh on Sh in (55) is a control for the 
continuous problem. 

Let us consider now the case of linear initial data given by a bi-grid algorithm, i.e. Sh is given by 
(56). Taking into account Propositions 6.2 and 6.3 we see that for initial data (U2,U^) £ Sh the Fourier 



representation ( 43 ) has the more particular form 



(89) 



(JV-l)/2 (JV-l)/2 
k=l k = l 



The basis function t/>f for the space Sh is given by 



(90) r h = & 



a , k W 4 (A a ' k ) 0tk W 5 (A a > N+1 - k ) a,N+i—k , W 5 (A^ N+1 - k ) W 4 (A a - N+1 - k ) , w+1 _ fe 



W±(A°>i 



W b {A a ' k ) 



-V h 



W 5 (A a > k ) W 4 (A°- N + 1 - k ) (fh 
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W 4 (A) 



A 



A - 10 



W(A) and W 5 {A) 



(60-A)(A-10)(A-12) 
(A 2 + 16A + 240) 2 



W(A) 



and W is defined by (39). Let us fix 1 < k < (N — l)/2 and consider the homogeneous problem (10) with 
initial data (U2,Uj^) = (s^/iXi' ,ip k ), s e {— 1,+1}, for which the solution takes the form 

/n-i \ tt /-t\ ^a.k/,\ a,k . -~-o.k/,\ o.k . -~~&.k/,\ a,N-\-l — k . ^o,fe/,\ o.N-\-\ — k 

(9!) Ufc(t) =u^ (t)<p h ' +u lo (t)(p h +u h ' i {t)<p h ' 

where the low frequency coefficients are 



s T7^fe exp(«aA h - (t-T)), (i) = 



V^ 4 (A°. fe ) ^ 2 



exp(±iAf (t-T)) 



and the high frequency ones are 

W 5 (A a ' JV+1 - fc ) 



^■a.k /,\ 
U M (*) 



and 



W 5 (A°. fc ) 
W 5 (A a > Ar+1 - fe ) ^(A '^ 1 -* 



± 



1 



. a,JV+l-fc 



exp(±iA? 



a,N+l-k 



(t-T)) 



W 5 {A 



a,k\ 



Wi{A° 



By considering the particular class of solutions given by (91) into (21), we see that the control Vhit) 
satisfies the identity 

T 



(92) 



Vh(t) exp(isX^' k i) dt 



{-If A" 



1 a,k 



where the error terms are 



£L = -(- 1 ) 



fe sin(fc7r/i) 



v h (t)(u°f(t)n°> k +u a h f(t)n 



~*~a,k i 



^a,N+l-k 



+ U M (*)« 



o,jV+l-fc 







and 

t-iA; ^o,k^o,k/r\\ i ^a,-/V+l — k^a,k / n \ . ^o,./V+l — k^o,k /r*\ ,k ,k / n \ ^a,./V-fl — k^-a,k /r\\ ^o,N+l — k^o,k /n\ 

E h,2 = yh,i u io (°)+yh,i u hi^) + v h ,i u h l { Q )-Vh,o u io,M-y h fi u hi,M-Vhfi u hl ,M- 

Passing to the limit as h — > in the left hand side and in the first term in the right hand side of ( 92 ) can 
be done as for the truncated space Sh in (55). Therefore, in order to prove that the weak limit v* satisfies 

h -t 0, i.e. 

as h -» 0. 



(87), it is enough to show that the error terms are small as h 0, i.e. 

(93) ->0 and \E k J 

From the fact that the L 2 (0,T)-norm of the discrete control Vh(t) is uniformly bounded as h 
Cauchy-Schwartz inequality and the explicit expressions of u°^ k (t), u c ^ k (t) and u^ k (t), we obtain 



0, the 



ith 



E 



k.l 

h,l 



,o,fc|2 



W 4 (A a ' k ) 2 /sin 2 (/br/i) sin^farft) 



W 4 (A°' k ) 



J^a,k 



and 



E, 



h'l = \n°' N+1 ~ k \ 2 



fe < 2 — |„o,JV+l-fe|2 



W 5 (A a ' N+1 - k ) 2 /sin 2 (fc7r/i) sin 2 (fc7r/i) 
W5(A a ^) V A^ A Q ' Ar+1 - fe 



W 4 (A^ Ar + 1 -' £ ) 2 T^ 5 (A Q ' Ar+1 -' £ ) 2 / S in 2 (fc7r/i) sin 2 (A:7r/i) 



W 4 (A 



o,N+l-k 



W b {A a ' k ) 



A a - 
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where W4 and W5 are as in (90 1. 

On the other hand, since the 
bounded as h — > 0, we see that 



norm of the initial data (Yi,Y°) to be controlled is uniformly 



K, 2 I < C(E K h 



k.l 



with 



< 2 = A, 



o.k \-^o,k 



77i/c,2 * a.N -\-l — k \ ^a,k 
E h2 = A h \ U hi 



~*~o,k 1 



w k,2 p fc,3\l/2 
2 _ W 4 (A a '^ 2 



(o)l 2 + K£(o)| 2 = 



Wi(A°' k 



A°>' 



1 



jya,k 

2 / ^a,N+l-k 
^a,k 



1 



and 



j^k,3 ao.AT+1- 

h h,2 = A h 



c (o)l 5 



l hi,t 



(o)l 2 - 



W 5 (A a ' k ) 

W4,(A a ' N+1 ~ k ) 2 W 5 {A°" N+1 - k ) 2 /A o,N+l-k 



Wi{A 



o.N+l-k\ 



W 5 (A a '< 



A a " 



1 



of the terms E 



k.i 



For a fixed k £ N, let us study the convergence as ft —r ^ U i u±ic . 
• From |n°< fc | 2 3TU(60) = 10 (W introduced in S), |W4(A Q ' fc )| 2 



, 1 < i < 3, 1 < j < 2: 

W 4 (A°< fe )| 2 - 
0. 

W^(0)| 2 



l^ 4 (0)| 2 = 0, 

|W 4 (60)| 2 = 24/25, sin 2 (fc7r/i)/A a ' fe -> 1 and sin 2 (/brft)/A°- fc -> 0, we see that E*'{ ^0as/i 
From |ri«.JV+i-fe|2 31 y( 10 ) = 0, |T^ 5 (A a ' iv+1 - fc ) | 2 -> | Ws(10)| 2 = 0, |W 5 (A Q < fc )| 2 
1/96 and sin 2 (fc7rft)/A a < jv+1 - fc -> 0, we see that £^> 2 -> as ft -> 0. 
Remark that | n °>W+i-fc|2 _^ 3^(12) = 6, | W 4 (A°' JV "+ 1 - fc )| 2 |WVl2)| 2 = 72, sin 2 (fc7rft)/A°- Ar+1 -' £ H 
0, but |W/ 4 (A Q < Ar+1 - fe )| 2 -> |VF 4 (10)| 2 = 00. Nevertheless, |W 4 (A a - Ar + 1 - fc )| 2 |W r 5 (A^ Ar+1 - fc )| 2 -> 0, so 
that at the end £^ as ft 0. 

Remark that A°- k 4 60 and A Q < fc -> 0. Observe that \W 4 (A a > k )\ 2 -> like |A a ' fc | 2 . This compensates 
the sing ularity of 1/A a ' fc so that E k ^\ -> as ft -> 0. 

| Wg ^ A o,JV+i-fe)|2 q since it invo i ves tne f actor ( 10 _ ^N+i-k^ ~ sin 6 (/brft/2). This cancels the 
singularity introduced by 1/A a,fc ~ 1/ sin 2 (fc7rft/2) and ensures that £^/ 2 ~~ * as ft — >• 0. 
|W 4 (A a -Ar+ 1 - fc )| 2 |T^ 5 (A a ^+ 1 -' £ )| 2 ^ q since it contains the factor (10 - A. a ' N+1 ~ k ) 2 ~ sin 4 (fc7rft/2) 
which compensates the singularity of 1/A a ' k ~ 1/ sin 2 (£;7rft/2) so that — > as ft — » 0. 



This concludes ( 93 1 and the fact that the weak limit u* of the sequence of discrete HUM controls obtained 
by minimizing the functional Jh over the bi-grid class Sh in ( 56 1 is a control in the continuous problem ([!]) . 

Step II - The weak limit v* is the normal derivative of a solution of the continuous adjoint 
problem ([3]). Consider (U°,U^) 6 Sh (which in what follows can be both the subspace in (551 or the one 
in (|56 1 ) to be the minimum of the functional Jh- Due to the uniform nature of the observability inequality 



(14 



£h(U2, Ui) is uniformly bounded, i.e. there exists a constant C > independent of ft such that 



(94) 



1 N 

-Y 

2 ^ 

k=l 



,k,0, 



,k,l, 



A 



o.k 



o,k.0, 



,k,l, 



< c. 



Due to property (57), the resonant mode in the solution of the adjoint problem (10) for initial data in the 
filtered space Sh in (55) or (56) vanishes, so that the Fourier representation of the total energy in the left 
hand side of (94) is valid for both filtered spaces Sh in (55) and (56). 

Remark however that the high frequency components vanish for data in the truncation subspace Sh in 



(55). On the other hand, for data in the bi-grid space Sh in (56), the relations between the optic and the 
acoustic modes and the high frequencies in the acoustic mode and the lowest ones described in Propositions 
[Ql and l6~3l hold. 



From (94), we deduce that 



;k,l 



(95) (A?'*fi 



, )k^{X k u ) fe , (u h ) k ^(u ) k , {-^) k , K )k 



as ft — >• 0, weakly in 
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Set u*' l (x) :— X^fcLi " ' ¥ k ( x )i observe that (w*'°, u*' 1 ) € V and denote by u*(x,t) the corresponding 
solution of ([3]). Firstly, we prove that 



(96) 



U N (t) 



"4(1, t) ash^O weakly in L 2 (Q,T). 



In fact, for arbitrary functions ip £ L 2 (0,T) and ip e £ Cq(0,T), we will prove the following estimate: 

T T 



(97) 



u N {t) 

h 



u* x (l,t)U(t)dt 



< 



r a u N (t) 
h 



u* x (l,t))ip(t)dt 



C\\iP-A\\LHo,T)+Ch k \\^\\ LHO:T)l 



where T a := T" is the projection on the acoustic branch defined by (72 1. In a similar way, we define the 
projection on the optic branch, T°. In order to prove (97 1, we decompose its right hand side as follows: 

T 



r a u N (t) 



-«:(m))v(*)#+ / ( 



r°u N (t) 
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Taking into account that ip e £ C k (0,T), by integration by parts, we have 
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From the Cauchy- Schwartz and the admissibility inequality, the bound (94) and the fact that A^' Z > 12 /h 2 
for all 1 < I < N, we deduce that 
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Once we get (97), we conclude (96) by using the following three ingredients: 
• the weak convergence (82) combined with the strong convergence n a ' 1 -> \/2, sm(lirh) / \ a < 1 -> 1 and 

T T 

^(i) exp(±i(i - T)X a h l ) dt -> J ip(t) exp{±i(t - T)X l ) dt as h -> 0, 



allowing us to pass to the limit as h — > in the sense of £ 2 in the sum in the right hand side of 



r a u N (t) 
h 



N 



^(i)*=EE A ^±^(-i) 



± 1=1 



sin(^7r/i) 



exp(±i(t-T)\ a h ' l )dt 



so that we can guarantee that the first term in the right hand side of (97) is small as h — > 0. 

the density of (0, T) in L 2 (0, T), allowing to choose e so that \\ip — '0c||l 2 (o,t) is arbitrarily small. 

an appropriate choice of the mesh size h according to e, so that h k \ \ipe | \l 2 (o,t) is arbitrarily small. 
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Let us check that v* — u*(l,t) in L 2 (0,T). Indeed, using (86) and then (96), we have 



T T T T 

v*(t)ip(t)dt = lim / v h {t)i>{t)dt= lim / ( - ^t^W) dt = I *£(MM*)<ft, 

h— >Q J h—tO J \ hi J 



for an arbitrary ip <E L 2 (0,T). Therefore, v* — v in L 2 (0,T) and also u*(l,t) = u x (l,t), which, jointly with 

gives that (u*'°, u*- 1 ) — (vPjU 1 ) in V and then (u ' ' )/. = (u ' )/- 
in £ 2 , i — 0, 1. This means that once we have identified that the weak limit of the discrete controls is the 
continuous HUM control, we have the T-convergence of the discrete minimizer to the continuous one. 

Step III - Strong convergence of the discrete controls. In order to prove that Vh converges 
strongly in L 2 (0, T) to v as h —> 0, it is enough to prove that 

T T 

(98) limy* \v h {t)\ 2 dt = J \v{t)\ 2 dt. 





Using as test function XJh(t) in (21) the minimizer Uh(t), we have that 



T T _ 
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where 



N ^a,k,l 



P h =2^ 2^2\~V^k =F fh J (*V u ^ ±u h ) exp(T*rA, i ' ). 



± fc=i - " iX h 



Let us remark that using the strong convergence of the acoustic part of the initial data to be controlled 



(82 ) and the boundedness of the energy of the minimizer of J~h, 

T T 



hmP^ = ^^V^-Ty fc ' )(^ fc ' ±^ 1 )exp( Ti TA' £ )= / \u x (l,t)\ 2 dt = / \v(t)\ 2 dt. 



± fc=i 



o o 



On the other hand, taking into account the uniform boundedness of the energy of the minimizer of Jy L1 we 
get 

N i — o, fe, 1 1 2 

\ P h\ 2 <C^ + ir/'V) as ^ 0, 



A 

fc=i 



which concludes (98) and the strong convergence of the optimal control. □ 



Remark 6. In [7], it was proved that, for initial data (y ,^/ 1 ) in the continuous control problem belonging 
to the more regular space Hq x L 2 (0, 1), the numerical controls obtained for the finite difference or the linear 
finite element semi- discretization of the wave equation |7]) converge to the continuous HUM controls with 
an error order ft, 2 / 3 . This is due to the fact that \\h(Q — £| ~ h 2 £ 3 , where \h(Q can be each one of the 
dispersion relations for the finite difference or finite element approximation of the wave equation. In |18j , we 
observed the fact that the acoustic dispersion relation A^(£) of the P% - finite element method approximates 
the continuous one £ with error order h £ 5 for all £ € [0, ir/h], so that the convergence error for the numerical 
controls obtained by the bi-grid algorithm in the quadratic approximation of the wave equation increases to 
/1 4 / 5 under the same regularity assumptions on the continuous initial data to be controlled. 

8. Final comments and open problems 
• All the results in this paper can be extended to finite element methods of arbitrary order k, with 



the additional difficulty that when computing the eigenvalues, the quadratic equation (30) has to 
be replaced by a fc-th order algebraic equation in A which is technically complicated to be solved 
explicitly. The same difficulty arises when passing to several space dimensions. 
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• The results in [5] providing a general method to obtain uniform observability results for time dis- 
cretizations of conservative system lead to the extension of our observability results for the P2 - space 
semi-discretization to fully discrete P2 conservative approximations of the wave equation. 

• The extension of the results in this paper to non-uniform meshes is a completely open problem. 

• The last open problem we propose is related to and concerns the dispersive properties of the 
Schrodinger equation approximated in space using the P2 - finite element method. Designing appro- 
priate bi-grid algorithms taking care of all the singularities of both group velocity and acceleration 
simultaneously is an open problem. 
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